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Using the color-singlet free energy Fi and total internal energy Ui obtained by Kaczmarek et al. 
for a static quark Q and an antiquark Q in quenched QCD, we study the binding energies and 
wave functions of heavy quarkonia in a quark-gluon plasma. By minimizing the grand potential in 
a simplified schematic model, we find that the proper color-singlet Q-Q potential can be obtained 
from the total internal energy Ui by subtracting the gluon internal energy contributions. We carry 
out this subtraction in the local energy-density approximation in which the gluon energy density can 
be related to the local gluon pressure by the quark-gluon plasma equation of state. We find in this 
approximation that the proper color-singlet Q-Q potential is approximately _Fi for T ^ Tc and it 
changes to -\-\Ui at high temperatures. In this potential model, the J/ip is weakly bound above 
the phase transition temperature Tc, and it dissociates spontaneously above 1.62rc, while Xc and 
, t/j' are unbound in the quark-gluon plasma. The bottomium states T, Xb and T' are bound in the 

quark-gluon plasma and they dissociate at 4.10rc, l.lSTc, and 1.38Tc respectively. For comparison, 
we evaluate the heavy quarkonium binding energies also in other models using the free energy _Fi or 
the total internal energy U\ as the Q-Q potential. The comparison shows that the model with the 
new Q-Q potential proposed in this manuscript gives dissociation temperatures that agree best with 
' those from spectral function analyses. We evaluate the cross section for a{g + J/t/) — > c -|- c) and its 

^ I inverse process, in order to determine the J/'tp dissociation width and the rate of J/ip production 

, by recombining c and c in the quark gluon plasma. 
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^ ' I. INTRODUCTION 

\ The stability of heavy quarkonia in the quark-gluon plasma is an interesting subject of current research in high- 
ip'' energy heavy-ion collisions as Matsui and Satz suggested that the suppression of Jlij) production can be used as a 
signature of the quark-gluon plasma Q. DeTar 0,13, Hansson, Lee, and Zahed 0, and Simonov H, IE 13 argued 
1^ ■ however that because the range of strong interaction is not likely to change drastically across the phase transition, 
low-lying mesons including J /if) may remain relatively narrow states and the suppression of J /tp is not a signature of 
the deconfinement phase transition Q . Whether or not J/'0 production will be suppressed depends on the screening 
between the heavy quark Q and the heavy antiquark Q when the quarkonium is placed in the quark-gluon plasma. 
The degree of screening is highly nonperturbative at temperatures near the phase transition temperature |g. The 
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' related question of the quarkonium stability must be examined in nonperturbative QCD using, for example, the lattice 
gauge theory. 

Recent investigations of masses and widths of h eavy quarkonia in quenched lattice QCD calculations were carried 
out by Asakawa et al. [ol. Il0l| and Petreczky et al. |lllll2l[T3l | using the spectral function analysis and the maximum 
entropy method. They found that the width of J/ijj remains relatively narrow up to 1.6 times the critical phase 
transition temperature Tc. Reconsidering the properties of the quark-gluon plasma also led Zahed and Shuryak to 
suggest that quark-gluon plasma at temperatures up to a few supports weakly bound meson states |l4l llSl \Wl Il7j . 
They have also estimated the binding energy of J/V' and found it to be stable up to 2.7 1}; . The possibility of weakly 
bound meson states in the quark-gluon plasma was suggested earlier by DeTar 0, Hatsuda, and Kunihiro |l8l |. 
Phenomenological discussions on medium modifications of charmonium in high-energy heavy-ion collision have been 
presented recently by Grandchamp, Rapp, and Brown |l9j |. Summaries of recent development in heavy quarkonium 
suppression and deconfinement have also been reported by Petreczky [S^ and Karsch [2l| . 

As the knowledge of the stability of J / ip has important implications on the fate of J/ip in the quark-gluon plasma, 
it is important to obtain an independent assessment on the binding of heavy quarkonia, in addition to those from 
previous analyses. The spectral function analyses of heavy quarkonia using gauge-invariant current-current correlators 
have been carried out in the quenched approximation. Within the quenched approximation, independent lattice gauge 
calculations have also been carried out using the correlation of Polyakov lines from which the free energy Fi and the 
total internal energy Ui can be calculated The two-body potential obtained from the lattice gauge theory 

can be used to study the dissociation of heavy quarkonia. It is of interest to ask whether, within the same quenched 
approximation, the spectral function analysis and the potential model analysis will lead to consistent results concerning 
the stability of quarkonia in the quark-gluon plasma. As we shall deal with lattice results from quenched QCD only. 
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the quark-gluon plasma wc shall consider consists of gluons. For convenience, we shall continue to refer to such a 
gluon medium from quenched QCD as a quark-gluon plasma. 

Besides checking the consistency of independent quenched lattice gauge calculations, we would like to use the 
potential model to examine many physical quantities of interest. If quarkonia are indeed stable in the plasma, it is 
useful to find out how strongly bound they are. Furthermore, J/V' can dissociate by collision with constituents of 
the plasma. In such a coUisional dissociation, the rate of dissociation depends on the cross section for the reaction 
g + J/4' c + c. We would like to evaluate this cross section as a function of temperature T, which can be obtained by 
using the bound state wave functions in the potential model. The knowledge of the dissociation cross section allows 
a determination of the coUisional dissociation width. 

In energetic heavy-ion collisions, many pairs of charm quarks and antiquarks may be produced in a single central 
collision [2^ . These charm quarks and antiquarks can recombine to form J/ tp in the quark-gluon plasma. We also 
wish to find out here the rate of producing J/ijj through such a reaction. The production rate depends on the cross 
section for the reaction c + c —f J/tp + g. The latter quantity can be obtained from the cross section for the inverse 
reaction g + J/ip ^ c + which we already intend to calculate. 

Previously, the effects of temperature on the stability of heavy quarkonium was studied by Digal et al.^^^^ and 
Wong [23,123 using the free energy and assuming that the effects of entropy are small. It was, however, pointed out 
by Zantow, Kaczmarek, Karsch and Petreczky |M l29ll30l| that the effects of entropy depend on the separation distance 
between c and c. They suggested that the total internal energy Ui , instead of the free energy Fi , may be used as the 
QQ potential for the calculation of heavy quarkonium bound states. As the theoretical basis for this suggestion has 
not been fully explained in the literature, we shall go into details to examine the theoretical questions on the proper 
potential for Q-Q states. We find that the proper Q-Q potential involves the Q and Q internal energy C/LU. We shall 

show that in the local energy-density approximation, = 3i^i/(3 + a) + aUi/{3 + a) where a = 3p/e is given by 
the quark-gluon plasma equation of state. 

When a heavy quarkonium is placed in a quark-gluon plasma, conventional description assumes that the medium 
effect is dominated by the effect of Debye screening 0, , which leads to a decrease in the attractive interaction 
between the heavy quark and antiquark. We would like to study the effects of antiscreening due to the deconfined 
gluons and the relationship between antiscreening and the area law of spatial Polyakov loops 0, 0j Isil Issl Is^ . We 
would like to show that because the Gauss law of QCD contains a non-linear term involving the gluon field, the 
gluon field induces color charges at the field points. These induced color charges act to antiscreen the interaction 
between the heavy quark and the antiquark. We shall show that the strength of the antiscreening effect increases 
with an increase in the gluon correlation length and is proportional quadratically to the magnitude of the gluon fields. 
The antiscreening effects due to deconfined gluons bring an additional degree of freedom to mediate the interaction 
between the quark and the antiquark. 

This paper is organized as follows. In Section II, we review the heavy quarkonium production mechanism and 
the thermalization of the quark-gluon medium in high-energy heavy-ion collisions. We examine the evidence for 
rapid thermalization as revealed by the elliptic flow and hydrodynamics. In Section III, we review the lattice gauge 
calculations for the interaction between a heavy quark and a heavy antiquark and the gauge dependence of the 
interaction in bound state problems. In Section IV, we show that the total internal energy Ui contains contributions 
from the internal energy of the Q-Q pair and the internal energy of the gluons. In Section V, we use an appropriate 
variational principle to obtain the equation of motion for the quarkonium single-particle states and find that the 
proper Q-Q potential involves only the Q-Q internal energy. In order to obtain the internal energy of the heavy 
quark pair, it is necessary to subtract the gluon internal energy from the total internal energy Ui. In Section VI we 
show how such a subtraction can be carried out in the local energy-density approximation, using the quark-gluon 
plasma equation of state and the First Law of Thermodynamics. In Section VII, we show how the color-singlet Fi 
and Ui obtained by Kaczmarek et al. in quenched QCD can be parametrized and the proper Q-Q potential can 
be obtained as a linear combination of Fi and Ui from the lattice gauge results. Using this heavy quark- antiquark 
potential, we calculate the eigenenergies and eigenfunctions for charmonia in the quark-gluon plasma as a function 
of temperature in Section VIII. The locations of the dissociation temperatures at which heavy quarkonia begin to be 
unbound are then determined. The heavy quarkonium dissociation temperatures are compared with those determined 
from spectral function analyses. In Section IX, we calculate the eigenenergies and eigenfunctions for bb bound states 
in the quark-gluon plasma as a function of temperature. We discuss the effects of antiscreening due to deconfined 
gluons in the quark-gluon plasma in Section X. In Section XI, we discuss how the J/ijj bound state wave function can 
be used to calculate the cross section for g -\- J/ip c-\-c after the J/ -0 absorbs an El gluon, using the formulation of 
gluon dissociation cross section presented previously . The dissociation cross sections and coUisional dissociation 
widths of J /ip in the quark-gluon plasma are then determined as a function of temperature in Section XII. In Section 
XIII we evaluate the cross section for the inverse process of c -I- c ^ J/i/' -f ,g using the cross section g + J ^ c + c 
obtained in Section XII. The rate of J /ip production by recombining c and c in a quark-gluon plasma is estimated. 
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Wc conclude our discussions in Section XIV. In Appendix A, we show that the integral of the gauge fields along a 
space-like Polyakov loop obeys an area law if the gauge fields are correlated. This result is used in Section X to 
explain the antiscreening effect. 

II. HEAVY QUARKONIA PRODUCTION AND THE THERMALIZATION OF THE MEDIUM 

We are interested in using a heavy quarkonium to probe the properties of the matter produced in central high- 
energy heavy-ion collisions. In the collider frame, the colliding nuclei have the shape of Lorentz-contracted disks. The 
collisions are known to be highly inelastic in which a large fraction of the incident collision energy is released after 
the collision. What is the rate of the relaxation of the initial configuration to thermal equilibrium? 

From the experimental viewpoint, recent RHIC experiments by the STAR PHENIX and PHOBOS ^ 
Collaborations reveal the presence of an elliptic collective flow in non-central Au-Au collisions at RHIC energies. 
The occurrence of such a flow indicates that the initial azimuthally symmetric momentum distribution of particles is 
deformed into an azimuthally asymmetric momentum distribution. The magnitude of the azimuthally asymmetry is 
sensitive to the time at which the free streaming of particles terminates and the dynamics of a thermally equilibrated 
system begins [s^ |3^ ^| . Too late a thermalization time will lead to a spatially more extended system with a lower 
pressure gradient and a smaller azimuthal asymmetry. The azimuthal asymmetry is also sensitive to the numbers of 
degrees of freedom in the equation of state. The magnitude of the asymmetry can be well explained in terms of a 
hydrodynamical model of the quark-gluon plasma by assuming thermalization at an initial time about 0.6 fm/c |38l |. 
We infer from the experimental elliptic flow data and its hydrodynamical description that the thermalization in the 
central region of a RHIC nucleus- nucleus collision is very rapid, as short as 0.6 fm/c after the collision. 

From theoretical viewpoints, it was first pointed out by Landau |4ll | that the initial configuration after a high-energy 
nuclear collision consists of matter at an extremely high energy density in a very thin disk. The great magnitude of 
the energy density means that the number density of quanta of matter is very large. Such a large number density in 
the thin disk of matter leads to a very small mean-free-path compared to its dimensions, leading to a rapid relaxation 
to thermal equilibrium. According to Landau, "in the course of time, the system expands, the property of the small 
mean-free path must be valid also for a significant part of the process of expansion and this part of the expansion 
process must have hydrodynamical character" . Landau hydrodynamics provides a reasonable description of the 
widths of the rapidity distribution for high-energy hadron-hadron and nucleus-nucleus collisions from ^/s = 3 GeV 
to RHIC collisions at 200 GeV 0| . Hydrodynamical description with a rapid thermal relaxation also provides a 
good description of the elliptic flow of matter after a Au-Au collision at RHIC, as indicated above. 

It is also useful to point out that a quanta in the non-equilibrium QCD matter interacts not only with other quanta 
(gluons and quarks) in two-body processes in terms of two-body collisions, but the quanta also interact with the fields 
generated by all other quanta. Because of the non-Abelian nature of the QCD interaction, the fields generated by 
other quanta are also sources of color fields and the quanta must in addition interact with the color fields generated 
by the fields of all other particles, in a highly non-linear manner (see Section X for another manifestation of the non- 
linear nature of the gauge field). Thus, a quanta interacts with other quanta not only by direct short-range two-body 
collisions, but also by highly non-local action-at-a-distance long-range interactions, through the fields generated by 
the fields of other quanta. There is thus an additional non-linear and long-range mechanism of thermalization in 
non-Abelian interactions which provides an extra push for rapid relaxation to thermal equilibrium. 

The rate of thermalization of a quark-gluon system after an ultra-relativistic heavy-ion collision is the subject of 
current theoretical research and has been discussed by Wong Molnar, and Gyulassy [isll. T he small mean- free- 
path has also been discussed by Shuryak Gyulassy and McLerran (see also Ref. |48l|'). The focus of the 
research is on trying to understand the phenomenologically fast rate of thermalization as indicated be the experimental 
elliptic flow evidence. For example, in the work of Molnar and Gyulassy |45| . it was necessary to shorten the parton 
mean-free-path by an order of magnitude in order to reproduce the magnitude of the elliptic flow. Similarly, in the 
work of Lin, Ko, and Pal [49l | in parton cascade, it was necessary to increase the parton-parton cross section by a 
large factor to describe the dynamics in a nucleus- nucleus collisions at RHIC. 

To use heavy quarkonia as a probe of the quark-gluon plasma, we need a knowledge of the heavy-quarkonium 
production mechanism. In a nucleus- nucleus collision at high energies, the partons of one nucleon and the partons 
of another nucleon can collide to produce occasionally a pair of heavy quark and antiquark. The time scale for the 
production is of the order of h/2mQ_ where mg is the mass of the heavy quark. It is of order 0.06 fm/c for a c-c 
pair and of order 0.02 fm/c for a b-b pair. As the initial partons carry varying fractions of the initial momenta of 
the colliding nucleons, the heavy quark pair will come out at different energies. Depending on the Feynman diagram 
of the production process, the produced Q-Q pair after the hard scattering process may be in a purely color-singlet 
quarkonia state or a coherent admixture of color singlet- and color-octet states [H^, l5lj |. The projection of different 
final states from a coherent admixture gives the probability amplitude for the occurrence of the final states. A color- 
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octet state will need to emit a soft gluon of energy -Egiuon to become subsequently a color-singlet state in an emission 
time of order fi./i?giuon- For the emission of soft gluon of a few hundred MeV, the time for the emission is of order 
0.5-1.0 fm/c. 

From the above considerations on the rapid thermalization of the quark-gluon plasma and the time for the production 
of heavy quarkonia, one envisages that by the time when the colliding matter is thermalized at about 0.6 fm/c, a 
large fraction of the quarkonia have already been formed, although in various energy states. The quark-gluon plasma 
is expected to have a life time of a few fm/c which is longer than the heavy quarkonium orbital period of order 
h/ (0.5 GcV). It is therefore meaningful to study the fate of a produced heavy quarkonia in a thermalized quark-gluon 
plasma at a finite temperature. The behavior of the heavy quarkonium system before the thermal equilibrium of the 
quark-gluon plasma and the interaction of a coherent Q-Q color admixture in the thermalized quark-gluon plasma 
are other topics which are beyond the scope of the investigation of the present manuscript. 

It should be pointed out that for a heavy quarkonium system in a quark-gluon plasma the thermalization of the 
quark-gluon medium does not necessarily imply the thermalization of the heavy quarkonium system. The former 
arises from the interaction among the light quarks and gluons, while the latter depends on the interaction between the 
heavy quarkonium and the constituents the quark-gluon plasma. Our evidence concerning the the rapid thermalization 
given above refers to the thermalization of the quark-gluon plasma and not necessarily the thermalization of the heavy 
quarkonium system in the quark-gluon plasma. If an isolated heavy quarkonium is placed in the thermalized quark- 
gluon plasma, the heavy quarkonium system is not in thermal equilibrium with the medium. It will interact with the 
medium as its density matrix will evolve with time. Given a sufhcient time that is longer than the heavy quarkonium 
thermalization time, the heavy quarkonium will also reach thermal equilibrium with its thermalized quark-gluon 
plasma. The thermalization status of a heavy quarkonium system can be inferred from the occupation number 
distribution of heavy quarkonium single-particle states. The occupation numbers in a thermalized heavy-quarkonium 
system will obey a Bose-Einstein distribution characterized by the temperature. In our present work, we shall study 
both a thermalized heavy quarkonium system and an isolated QQ bound state in the quark-gluon plasma. 



III. LATTICE GAUGE CALCULATIONS 



In a quark-gluon plasma, a quarkonium is actually a heavy quark and a heavy-antiquark each surrounded by a cloud 
of gluons and quarks. In the quenched approximation in which there are no dynamical quarks, the cloud surrounding 
the heavy quark and antiquark is approximated to consist of gluons only. As gluons are involved, the quark-antiquark 
system will be in different color states at different instances. We shall be interested in those systems in which the heavy 
quark and antiquark exist in the color-singlet state. Only in the color-singlet state will be the effective interaction 
between a quark (plus its cloud) and an antiquark (plus its cloud) be attractive. Such a color-singlet system can 
further absorb a gluon and become a color-octet system and we shall also study the cross section for such a process. 

The interaction between a heavy quark and a heavy antiquark in the color-singlet state was studied by Kaczmarek, 
Karsch, Petreczky, and Zantow |22|. They calculated {trL{r /2)L^ {—v /2)) in the quenched approximation and they 
obtained the color-singlet free energy Fi{v,T) from 

(iri(r/2)Lt(-r/2)) = e-^i("-'^)/'=^. (1) 

Here <rL(r/2)Lt(— r/2) is the trace of the product of two Polyakov lines at r/2 and — r/2. The quark and the antiquark 
lines do not, in general, form a close loop. As a gauge transformation introduces phase factors at the beginning and the 
end of an open Polyakov line, {trL{v /2)L'^ {—t /2)) is not gauge invariant under a gauge transformation. Calculations 
have been carried out in the Coulomb gauge which is the proper gauge to study bound states. 

It should be noted that while the interaction between the quark and the antiquark is gauge dependent, the bound 
state energies are physical quantities and they do not depend on the gauge. As we explain below, a judicial choice 
of the Coulomb gauge in the bound state calculation will help in avoiding spurious next-to-leading contributions and 
singularities, which can be removed in other gauges only by additional laborious effort [S^, IH, Im ]. 

To study the bound states of a heavy quarkonium, we need a bound-state equation, such as the Bethe-Salpeter 
equation, and the interaction kernel in the equation. The non-relativistic approximation of the Bethe-Salpeter equation 
leads to the usual Schrodinger equation with the gauge-boson-exchange interaction [ssl Is^ . It is necessary to choose 
a gauge to specify the gauge-boson-exchange interaction. We can consider the case of QED from which we can get 
a good insight on the gauge dependence. For the static non-relativistic problem, the natural choice in the gauge- 
boson-exchange potential is the Coulomb gauge, in which the 1/q^ behavior is found in single Coulomb photon 
exchange. The binding energy, which is of order a^, has corrections only in the a* order. It gives the correct Breit 
equation with the proper spin properties when we expand the interaction to the next order. Graphs with the cross 
two-Coulomb-exchange diagrams vanish in the static limit, and uncrossed multiple Coulomb exchanges are strictly 
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iterations of the potential |52 |. In any other gauge, the zero-zero component of the photon propagator has some 
residual non-instantaneous contributions. A large number of Bethe-Salpeter kernels need to be included to eliminate 
the spurious contributions in the next-order of a'^ and a^/\na corrections [s^. l54j . Therefore, in their work on the 
static potential in QCD, Appelquist, Dine, and Muzinich suggested that the gauge freedom can be used to eliminate 
spurious long-range forces at the outset. They found that the Coulomb gauge continues to be useful in the static 
potential in QCD. The dynamics is now considerably complicated but spurious contributions are still eliminated |52l |. 

Based on the above, it is therefore important to recognize, as in QED, that there is a gauge-dependence in the 
two-body bound-state potential in the Bethe-Salpeter equation but it is most appropriate to solve the bound state 
problem using two-body potentials obtained in the Coulomb gauge, as was obtained by Kaczmarck et al |22| . 

IV. HEAVY QUARKONIUM STATES IN A THERMALIZED QUARK-GLUON PLASMA 

The state of a heavy quarkonium in a quark-gluon plasma can be described by a density matrix. The set of 
single-particle states for this density matrix can be chosen such that they can be represented well by quarkonium 
states in a Q-Q potential, and the residual interaction between the gluons with the quarkonium can be treated as a 
perturbation. In this single-particle basis, the heavy quarkonium density matrix can be approximated to contain only 
diagonal matrix elements representing the probabilities for the occupation of different single-particle states. What is 
the Q-Q potential that enters into the Schrodinger equation for these quarkonium single-particle states? 

The Q-Q potential in perturbative QCD has been studied by Petreczky We would like to examine here the 
Q-Q potential in non-perturbative lattice QCD calculations. In previous analysis, the Q-Q potential was taken to be 
the free energy Fi for a pair of correlated Polyakov lines, assuming that the effects of entropy are small p^L l27i l28l | . 
It was, however, pointed out by Zantow, Kaczmarek, Karsch and Petreczky [l,!!^!!^ that the effects of entropy are 
large and the total internal energy J7i, may be used as the Q-Q potential. As the theoretical basis for these suggestions 
has not been fully discussed in the literature, we shall go into details on the proper description of the Q-Q potential 
and single-particle states. 

We start first by studying the auxiliary problem of a static color-singlet Q-Q pair at a separation r in a thermalized 
quark-gluon plasma. In the quenched approximation, the color-singlet free energy i^i(r,T) for such a static pair, can 
be written from Eq. ^ explicitly in the Euclidean time r = as |3lj 

e-f'^-^'^'^^ = Z{v,T) I Z,,{T), (2) 
Z{v,T) = j[dA]WQQ{A,T,v), (3) 
Zo(T) = j [dA]Wo{A,T), (4) 



Wqq{A,T,v) = tr <; Pexp{ / .g(r)drAo(-, r)} exp{ / .g(r)dTylo(--, r)} ) Wo{A,T), (5) 



^ J fj 



Vro(A,T) =exp|-i^''dT j d^xF^.F^-j, 



(6) 



where (3 = 1/kT is the inverse temperature, Z{r, T) is the partition function when a color-singlet Q and Q separated 
by a distance r is placed in the gluon medium, Zq{T) is the partition function in the absence of Q and Q. The 
operator P is the path-order operator, which is the time-order operator T for exp{J^ g{T)dTAo{r/2,T)} and is the 
reverse time-order operator for exp{/^ 5(r)(irylo(— r/2, r)}. The free energy Fi(r,T) with the Q-Q pair is measured 
relative to the free energy Fo{T) without the Q-Q pair. We rewrite Eq. jSJ as 

Zo(T) = J [dA]e^^^^'^-^hrlpeMj\{r)dTAo{^,T)}exp{j\{T)dTAo{^^^ J d'xF^^F^A .(7) 
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Taking the derivative of this equation with respect to (3, we obtain for the derivative of the left-hand side 

a{LHS} _ dZoiT) 



= I [dA] |-i I d'xF^,F^'-'^WoiA,T), 



(8) 



= j [dA]eP^^^-^) i^F,{v,T) + P^-^^^-\j d'xF,^F^-^WQQ{A,T,v) 

(9) 



d(3 d(3 
and for the derivative of the right-hand side, we get 

a{RHS} 

dp 

+tr|g(r)(Ao(^,r)-Ao(-^,r))Pexp{ 5(r)drAo(^,T)}exp{ g(r)dr (- ^, r)} [> W^o( A 

We equate the above Eq. JHJl to Eq. ©. Using e^-'^i^'''^) ~ Zo{T)/Z{r,T) and dividing the resultant equation by 
Zo(T), we obtain then the proper thermodynamic equality relating Fi, Si, and Ui, for a system with a color-singlet 
Q and Q separated r at temperature T, 

Fi{r,T) + TSi{r,T) = Uiir,T), (10) 

where S*! (r, T) ~ —dFi (r, T) / cJT is the color-singlet entropy with the Q-Q pair and is measured relative to the entropy 
So{T) = —dFo{T)/dT without the Q-Q pair, and Ui{r,T) is the total color-singlet internal energy given explicitly by 

Ui{r,T) = U^'Ur,T) + U^'\v,T) - Ugo{T), (11) 



U,oiT) = J [dA] d'xF^^F^''^ WoiA,T) j J [dA]WoiA,T), 

C7(i)(r,r) = J [dA] |i J d'xF.^F'^-'^ Wqq{AT,v) / J [dA]WQQ{AT,v) 



(12) 



(13) 



and 



U^Ql{r,T) = J [dA]tr}^g{T)[AoC-,T)-Ao{-l,T)]Pcxp{j\{T)drAo{'^,r)}c^^^ 

\dA]tr\pcxp{ [\{T)dTAo{^,T)}cxp{ ^ g{T)dTAo{-^,T)}\ Wo{A,T), (14) 



I Jo ^ J/3 J 

We may attempt to give names to various mathematical expressions. In Euclidean time, the quantity Ff_ii,F^^ / A is 
equal to (E^ + B^)/2, the gluon energy density [sills^ . The quantity Ugo(T) is the expectation value of / d"^x{E^ + 
B^)/2 with the weight function Wo{A,T), and it corresponds to the gluon internal energy in the absence of the 
heavy quark pair. It is independent of the separation r between Q and Q. In contrast. Ug^\r,T) is the expectation 
value of / d'^x(E'^ + B-^)/2 with the weight function WQQ{A,T,r), and it corresponds to the gluon internal energy 

in the presence of the heavy quark pair. Consequently, Ug^\r,T) depends on the separation r between Q and Q. 
The difference between the total internal energy Ui and gluon internal energy difference Ug^\r,T) — Ugo{T) is the 
quantity Uqq{t,T), the internal energy of the heavy-quark pair, including the interaction between Q and Q as well 
as the interaction between Q with gluons and Q with gluons. 

Eqs. (|10|) - (|14|l show that the total internal energy U^^^ (r, T) contains the gluon internal energy difference U^^^ (r, T) — 
UgoiT). In order to obtain the r-dependence of the internal energy of the heavy quark pair uj^Q{r, T), it is necessary 

to subtract the gluon internal energy difference C/f ^ (r, T) - C/go(T) fr om the total internal energy C/i(r, T). In Section 
VI we shall show a method to carry out such a subtraction in the local energy-density approximation. 
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V. EQUATION OF MOTION FOR Q-Q SINGLE PARTICLE STATES 

Lattice gauge calculations provide information on the free energy Fi and the total internal energy Ui for a static 
color-singlet Q and Q separated by a distance r. Quantities for Q and Q in the color-octet state can be similarly 
obtained. For simplicity, we shall limit our consideration to a system of color-singlet QQ states. The generalization 
to a system color-octet states can be easily carried out. 

We would like to use an appropriate variational principle to obtain the equation of motion for the color-singlct 
quarkonium states in a quark-gluon plasma. For such a purpose, we consider a schematic toy model that retains 
the relevant features of the system. In the quenched approximation without dynamical light quarks, the quark-gluon 
plasma consists of gluons only. The quark Q and antiquark Q are in dynamical motion in different single-particle 
states. For simplicity, we presume that the color degree of freedom has been integrated out. The dynamical degrees of 
freedom in our schematic model are then the QQ and gluon states, ipi(r) and 4>i(v), which can be boimd or imboimd, 
the corresponding QQ and gluon state occupation numbers, ni{QQ) and nj{g), and the total number of gluons Ng. In 
the schematic toy model, we represent the Q-Q, Q-g, Q-g, and g-g interactions when Q and Q belong to a color-singlet 
state by Vqq, Vqg, Vg^, and Vgg respectively. 

The equilibrium of the system at a constant temperature T is characterized by minimizing the grand potential A 
appropriate for the color-singlet Q and Q in dynamical motion in the quark-gluon plasma given by 

A = Ti- ii{QQ)NQQ - ii{g)Ng =Ui- TSi - KQQ)Nqq - K9)Ng, (15) 

where J-i is free energy, Ui the total internal energy, and Si the total entropy for Q and Q in dynamical motion in 
the quark-gluon plasma. Only in their static limits when the heavy quark and antiquark are held spatially fixed can 
J-'i, lAi and Si be equal to their corresponding static thcrmodynamical quantities Fi, Ui, and Si. The quantities 
fi{QQ) and fi{g) are the chemical potentials of QQ and gluons respectively, and Nqq and Ng are the numbers of 
QQ and gluons respectively. Strictly speaking the number of gluons at thermal equilibrium depends on the length 
scale (the Q^ value of the measuring probe) under consideration. In our schematic model, we can fix the length scale 
appropriate for Q-Q bound states in the quark-gluon plasma, and in that length scale the number of gluons, for the 
fixed spatial volume under consideration at thermal equilibrium at T, can be determined. We shall also ignore the 
annihilation of Q-Q into light quarks or photons and the corresponding inverse production so that the number of QQ 
pair can be considered fixed also. In the above equation, we need to add the Lagrange multipliers Xi{QQ) and Xj{g) 
to the grand potential in order to constrain the normalization of the wave functions. 

To carry out the minimization of the grand potential to obtain the equation of motion of the single-particle states, 
we can follow Bonche, Levit, and Vautherin [S^, III, [5^ |60, |6ll and write down the grand potential A exphcitly. The 
internal energy Ui in Eq. ^ is U^^^ +U^g^^ - Ugo, the sum of the internal energy of the heavy quark pair and the 
gluons, relative to the gluon internal energy when the heavy quark pair is not present. In terms of the wave function 
and occupation number degrees of freedom, the grand potential A can be written explicitly as 

A = 5]n.(gg) / drV'l(r)[|^ + yQQ(r)]V',;(r)+^n,(QQ)n,(g)(^,</.,|yQg + T/Qj^,^,) 

+ / '^^''l^]i^')\/pl+^ls + J2^i'^9)nki9){<Pj4'k\Vgg\(t)j(t)k + (t)k<Pj)/2 ~ Ugo 

+ T ^ [n,, {QQ) In n, {QQ) - {l + n,{QQ)} \n {l + n,{QQ)}] 

i 

+ T ^ [n, {g) In n, {g) - {1 + n, {g)} In {1 + {g)}] 
j 

- KQQ) E - /^(ff) E '^J- (-9) - E (V'. 1^^) - E (a) (16) 

i j i j 

Here p and /ircd are the relative momentum and the reduced mass of Q-Q, Pg and rrics are the momentum and the 
effective mass of the gluon, and the dependence of various quantities on the temperature is made implicit. In this 
expression for the grand potential, the first two terms give I^qq, the third and fourth terms give h(g^\ and the sixth 
and seventh terms give the entropy TSi. The matrix element {4'i4>j\^Qg + ^Qgl'^i'Pj) 



{ij.cj,,\VQg + VQgH'.,cj,,) = I drdr'^l{r)<j>l{r')[VQg{r' + ^) + ^^^(r' - ^)]V',(r)</>,(r'), 



(17) 
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and the matrix element {4>j4'k\Vgg\4>j4>k + 4>k<t'j) is 

(</>,0fc|ygJ0,(/<, + 0,0,) = / dr'dr''</>t(rO0l(r'')K,3(r'-r'')[0,(r')0fe(r'') + </'fc(r')'/>,(r'')]. (18) 



When we carry out the minimization of the grand potential with respeet to the dynamical degrees of freedom, we 
obtain five equations (Eqs. (24)- (28) below). By minimizing with respect to ipliy), we obtain the equation of motion 
for ?/',;(r) of the QQ system as 

'''P' +yi(r,r)-e:;(QQ)U.(r)=0, (19) 



where the color-singlet single-particle potential Vi{y,T) is 

yi(r,T) = FQQ(r) + ^,^,(5) / dr>](r')[yQ,(r' + r/2) + ^^^(r' - r/2)], (20) 

and e[{QQ) = Xi{QQ)/ni{QQ). We note that in the above equation, the quantity 

5]n,(g)^](r')0,(r')=P,(r') (21) 
j 

is the gluon density Pg{r'). The color-singlet potential Vi{r,T) is 

V^ir,T)^VQQir,T)+ J dr'p^lr', r)[FQ,(r' + ^) + y^^Cr' - |)], (22) 

which represents the internal energy of a static color-singlet Q and Q at a separation r in the quark-gluon plasma. 
It has the same physical meaning as internal energy uj^)^{r, T) of a static color-singlet Q and Q at a. separation r in 

the lattice gauge calculations. It is therefore reasonable to identify the QQ internal energy U^^Uv.T) of the lattice 

____ WW 

gauge calculations as the single-particle potential Vi(r,T) in the Schrodinger equation (|19fl . and write it as 



+ U)^l^{v,T)~e[{QQ)\i.,{v)^Q. (23) 



It is simplest to re-calibrate the energy (and similarly, the chemical potential) by ^i{QQ) = t[{QQ) — t^QQ(r)r-»oo 
and re-write the above equation as 

1^ + u'^qI{v,T) ~ f/W (|r| cx),T) - £,(QQ)| ^'.(r) = 0. (24) 

The above Schrodinger equation involving U^^q{v^T) is the proper equation of motion for quarkonium single-particle 
states. 

From the above considerations, the Q-Q potential in the equation of motion of quarkonium single-particle states 
is uj^)^{r,T). Lattice calculations so far provide information only on the free energy Fi(r, T) and the total internal 

energy J7i(r, T). It will be of great interest in future lattice work to evaluate J7g^(r, T) so that it can be used as the 

proper Q-Q potential in quarkonium studies. In the next section, we shall present a method by which J7^'^(r,T) can 
be approximately evaluated by using the quark-gluon plasma equation of state. 

The minimization of the grand potential with respect to the gluon wave function 0^ (r) gives the equation of motion 
for the gluon states. 



/p^+^c'ff + E ^» / '^'•^I i^Qs + '■Z^) + (r' - r/2)] V. (r) - (5) | (r') 
Y^n.ig) [ dr'Vt(r")F,,(r"-r')[0.(r")<^,(r')+</'fe(r')</',(r")] 



T.^x{9)nk{g) J dr>t (i.^)^t (^//) g^g.(jj^ [0,(r")0,(rO + Mr')Mr")]n,{9)Mr') = 0, (25) 



A,fe 
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where we have taken the density-dependent interaction as a deUa-function in r" — r' as in the Skyrme interaction 
|63l l64| ■ The minimization of the grand potential with respect to ni{QQ) and nj{g) yields 

"'^Q'^) = e[^.(QQKM(QQ)]/T_i ' (26) 

and 

^^•(g)= e[-.to)-Mg)]/T,i ' (27) 

which are the well-known Bose-Einstein distributions. 

The variation of the grand potential A of Eq. (15) with respect to Ng gives 

dA/dNg = /i(g). 

The requirement that thermal equilibrium is reached when the variation of A with respect to Ng is a minimum leads 

to m 

Kg) - 0. (28) 

The number of gluons at level j, nj{g), can then be obtained from Eq. (27) with //(g) ~ 0, and the total number of 
gluons Ng is given by 

^9 = E"^C9). (29) 

j 

The above considerations give the set of equations H24|l . H25|) . (|26|l . H27|) . and (28) for a system in which both the 
quark-gluon plasma and the QQ are in thermal equilibrium. If a quarkonium is placed in a thermalized quark-gluon 
plasma for a period longer than the time needed for it to thermalize, the heavy quarkonium will reach thermal 
equilibrium and the QQ system will be described by the above set of single-particle states with the Bose-Einstein 
distribution of occupation numbers. 

Other cases of our interest are those in which the quark-gluon plasma has reached thermal equilibrium while the 
QQ system, which arises from the independent mechanism of nucleon-nucleon hard scattering, may not have. For 
such a case, the equations of motion for the gluons states and gluon occupation numbers, Eqs. H25|l and 127|1 . remain 
valid, except that in Eq. H25|l the QQ occupation numbers ni{QQ) will no longer obey the Bose-Einstein distribution 
(|27|l . The equation of motion for the QQ single-particle states, Eq. (|24|l . remains valid, as it depends on the gluon 
density. Because the heavy-quark pair is a rare occurrence and gluons are much greater in number, we expect that 
the gluon density and wave functions obtained in Eq. (|25|) is insensitive to the status of QQ thermalization, and they 
depend mainly on the thermalization of the gluons themselves. Therefore, single-particle states of the QQ system, 
'4'i{v), which depend on the gluon density as given by Eq. H23|l . are relatively insensitive to the thermalization status 
of QQ, and they depend mainly on the thermalization status of the gluons. 

In view of the above, the set of QQ single-particles states il)i{v) of Eq. H24|) can be used to examine the states of a 
QQ system, whether the QQ has reached thermal equilibrium or not. For example, one can introduce a bound QQ 
state ^\ into a thermalized quark-gluon plasma of temperature T . Such a system is described by a density matrix 
with an initial occupation number ni(QQ) = Six at t = 0. The collision of the QQ with gluons in the medium will 
lead to the evolution of the occupation numbers of the QQ states as a function of time, leading eventually to the 
Bose-Einstein distribution characterized by the temperature of the medium. 

The QQ system can exist in color-singlet and color-octet states. The consideration we have given can be generalized 
to color-octet states. The Schrodinger equation for the different single-particle color states a will depend on J7^"^(r, T). 

For the same reason that the QQ single-particle states depend on the thermalization status of the gluons and are 
insensitive to the thermalization status of QQ, these single-particle states can be used to examine a QQ system in 
thermalized gluon matter. When the QQ also reaches thermal equilibrium, the set of states with occupation numbers 
given by equations will include both color-singlet and color-octet states. 

VI. COLOR-SINGLET Q-Q POTENTIAL 

From the variational principle of minimizing the grand potential, we find that the Schrodinger equation (|24|) contains 
the Q-Q internal energy U^q as the quarkonium Q-Q potential. Lattice calculations so far provide information only 
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on the free energy Fi(r, T) and the total internal energy Ui{r) but not t/lU. While we await accurate lattice gauge 

results for Uqq, we can obtain an approximate Uqq from Fi and C/i by using the equation of state of the quark-gluon 

plasma and the First Law of Thermodynamics. 

In the quenched approximation, the quark-gluon plasma consists of giuons only. We shall use the following strategy 

to obtain From lattice gauge calculations, we can calculate Fi — Ui which is equal to the entropy content TSi 

WW _ 

of the whole system. As the heavy quark Q and antiquark Q are held fixed in the lattice calculation, the entropy 
TSi comes entirely from the deconfined giuons, and it contains no contribution from the heavy quark QQ pair. From 
the entropy content of the gluon medium T^i, we can deduce approximately the gluon internal energy Ug^^ by using 
the gluon medium equation of state and the First Law of Thermodynamics. Then, the knowledge of the approximate 
gluon internal energy ujj^^ and the total internal energy Ui gives a color-singlet potential uj^Q in terms of Fi and Ui. 

Following such a strategy, we consider a deconfined gluon medium at temperature T with a static color-singlet Q 
and Q separated by a distance r. We focus our attention on a volume element dV of giuons at spatial position x in 
which the gluon internal energy element is dUg^^ and the gluon entropy element is dSi. At this spatial location, a gluon 
experiences an external potential V(x) exerted by the color-singlet Q and Q, other constituents, and induced charges. 
Hydrostatic equilibrium is reached when the local gluon pressure p^^^ (x) at the volume element counterbalances the 
external forces. We can begin in the non-relativistic description in which the total external force acting on a unit 
volume at x is pg(x)VxV(x) and hydrostatic equilibrium is reached when |6!tI | 

Vxp(i)(x) + pg(x)VxV(x) = 0. (30) 

In the relativistic case, one can describe the external interaction V(x) as part of the g^'^ metric, g'^'^ = 1 + 
2/9g(x)V(x)/e(x), and hydrostatic equilibrium is reached when |6^l6^ 

Vxp(iHx) + {p^^Hx) + 6(i)(x)}Vx\/l + 2p,(x)V(x)/6W(x) = 0. (31) 

Thus, the external forces determine the gluon pressure p^-^"^ (x) that is required to maintain hydrostatic equilibrium. 
At thermal equilibrium, the local gluon pressure is related to the local gluon entropy and the local gluon energy 
density. The external force therefore leads to a change of the gluon local entropy and the local gluon energy density. 

According to this picture, the change in entropy is zero when the heavy quark sits on top of the heavy antiquark 
at r = 0, as the color charges will neutralize for a color-singlet QQ. As Q separates from Q, the local gluon pressure, 
energy density, and entropy will increase in the vicinity of Q and Q to counterbalance the forces due to Q and Q. 
The increase in entropy will reach a constant value at large separations when the Q and Q each independently causes 
a modification of the gluon distribution in its vicinity. 

Under the temperature T and pressure p^^^(x), the gluon internal energy element dllg^"^ and the gluon entropy 
element dSi at x are related by the First Law of Thermodynamics, 

dUjj^'> = TdSi - {x)dV, (32) 

where the superscript (1) denotes that the heavy quark QQ pair is in a color-singlet state. The local gluon internal 
energy density eg^^(x) is therefore related to the local gluon entropy density dSi/dV and local gluon pressure p^^\yL) 
by 

eW(x) = ^(x)=r|i(x)-pW(x). (33) 

The equation of state of a homogeneous quark-gluon plasma (gluon medium only in the case of the quenched approx- 
imation) has been obtained in previous lattice calculations j68lj . It can be represented by expressing the ratio p/(e/3) 
as a function a{T) of temperature, 

a{T) = ^^. (34) 

We shall make the local energy-density approximation in which the local gluon energy density e'p (pi) and the local 
gluon pressure p^^)(x) under the temperature T obey the equation of state for the (homogeneous) bulk quark-gluon 
plasma so that 



_pW(x)_ 

4'^(x)/3 



a(T), (35) 
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as in the usual hydrodynamical description of the quark-gluon plasma. We then have 

dUg^^ . 3 TdS\ 

^ 3 + a(r) dV (''^ 

When we integrate the local gluon energy density over the whole volume of x under consideration, we obtain 

In a lattice gauge calculation, the large degrees of freedom of the gauge field and the maintenance of a thermal 
bath of constant temperature T make it reasonable to assume local thermal equilibrium and use the equation of 
state of homogeneous bulk matter to relate the local energy density to the local pressure in the local energy-density 
approximation. It will be of great interest to check the above equation H37|l in future lattice gauge calculations in 
order to test the validity of the local energy-density approximation. 

For the case in the absence of the Q-Q pair, the gluon internal energy Ugo{T) and the gluon entropy So{T) are 
similarly related by 

UgoiT) = ^^l^j.f So{T). (38) 

We therefore have 

C7W(r,r) - Ugo{T) = _A_T{5i(r,T) - SoiT)} . (39) 

In the above discussions from Eq. (32) to Eq. (39), we have followed the standard practice to use the convention of 
"absolute" units in which the entropy Si (absolute) and 5*0 (absolute) are measured relative to zero entropy. On the 
other hand, in lattice calculations and in Section IV, we have used the "lattice gauge calculation" convention in which 
the entropy 5*1 (lattice) and the free energy Fi (lattice) with the Q-Q pair are measured relative respectively to the 
entropy 5*0 (absolute) and free energy Fo (absolute) without the Q-Q pair (i.e. 5*1 (lattice) = 5*1 (absolute) — S'o (absolute) 
and Fi (lattice) = Fi (absolute) — Fq (absolute)). Henceforth, we shall switch back to the "lattice gauge calculation" 
convention in which the entropy 5*1 (r, T) and the free energy Fi(r, T) with the the Q-Q pair are measured relative to 
those in the absence of the Q-Q pair. In this lattice gauge convention, the above equation H39|l can be re-written as 

C/(i) (r , T) - Ugo {T) = ^-^^TS, (r , T) . (40) 

Using the relation TSi = Ui — Fi of Eq. (10) in lattice gauge theory, we can express the above Ug^\r, T) — Ugo{T) 
in terms of Fi(r, T) and Ui{r, T), 

C/W(r,r) - UgoiT) = _A_{C/i(r,T) -Fi(r,r)}. (41) 
Substituting this equation in Eq. (|11|) . we obtain 

The Q-Q potential which appears in the Schrodinger equation for quarkonium bound states, Eq. (|24|l . is then 
U^Qliv.T) - U^^lir cx),T) = fpiT) {F,{r,T) - Fi(r ^ oo,r)} + fuiT) {Ui{v,T) - C/i(r ^ (X3,T)}, (43) 



where 



iFiT) = (44) 
3 + a{l ) 



fuiT) 



3 + a{T) ' 



(45) 
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FIG. 1: (a) The energy density and pressure of a SU(3) gauge theory as a function of the temperature obtained by Boyd et al. 
Ig^ . (b) The ratio of 3p/e as a function of T/Tc. (3) The weight fp of _Fi and the weight fu of Ui that comprise the Q-Q 
potential according to Eq. ll43 



and fpiT) + friT) ~ 1. We shall use such a relation to obtain an approximate Q-Q potential from Fi and Ui. 

The above Eq. H43|) has been obtained for a color-singlet QQ pair in the quenched approximation. It can be easily 
generalized to the case of the unquenched full QCD with dynamical light quarks. In that case, the function a{T) 
corresponds to p/(e/3) appropriate for the equation of state of the quark-gluon plasma under consideration, and 
Ugo{T) of Eq. (|42|l becomes Uqgp{T), the internal energy of the quark-gluon plasma in the absence of the QQ pair, 
and Eq. (|43|) remains unchanged. 

We show e/T^ and p/T^ obtained in quenched QCD by Boyd et al. [b^I as a function of T/Tc in Fig. 1(a). In Fig. 
1(b) the function a{T), defined as the ratio 3p/e, is plotted as a function of T/Tc- The free energy fraction friT) 
and the internal energy fraction fu{T) calculated using this ratio of a{T) are shown in Fig. 1(c). One finds that at 
temperatures close to Tc, fp is close to unity, and the Q-Q potential is close to the free energy -F'i(r). As temperature 
increases, the Fi fraction decreases but approaching fp ^ 0.75 at very high temperatures. The Ui fraction is nearly 
zero at temperatures near Tc and it increases monotonically as a function of temperature, reaching a value of 0.25 at 
very high temperatures. 

It is of interest to discuss the conditions under which the application of the static potential uj^)=.{r,T) and its 
representation in terms of Fi and Ui as given in Eq. (|43|) can be reasonable concepts. For the static potential and 
the quark-gluon plasma equation of state to be applicable, it is necessary that the time for the quark-gluon plasma 
to reach thermal and hydrostatic equilibrium is short compared with the time for the periodic motion of the Q and 
Q. The time for the quark-gluon plasma to reach thermal equilibrium is of the order 0.6 fm/c, as one may infer 
from discussions in Section II. The orbiting time for a heavy quarkonium is of order 2r,-ms/v where r^ms is the root- 
mcan-square radius of the heavy quarkonium system and v is the relative velocity which is at most of order 0.5 for 
heavy quarks. The spatial scale of the heavy quarkonia in the quark gluon plasma is quite large. As we shall see in 
Section VIII and IX rmis = 0.88 fm at T/Tc = 1-13 for J/tp and it increases to rrms = 5.3 fm at T/Tc = 1.65. For 
Tf, ''rms = 0.25 fm at T/Tc = 1-13 and it increases to rmis = 0.59 fm at T/Tc = 2.6. The large spatial scales arises 
because these heavy quarkonium states are basically only weakly bound. The orbiting time for heavy quarkonia in 
the quark-gluon plasma is much greater than the quark-gluon plasma thermalization time. It is therefore reasonable 
to use a static Q-Q potential and the equation of state of the quark-gluon plasma to study heavy quarkonium in the 
quark-gluon plasma. 
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VII. SIMPLE PARAMETERIZATIONS OF Ui AND F-^ 

Kaczmarek et al. obtained the color-singlet free energy Fi(r,T) and internal potential Ui{r,T) in quenched 
QCD as a function of r = |r| and T. The radial dependences of Fi(r,T) and Ui{r,T) in quenched QCD can be 
adequately represented by a screened Coulomb potential with a screening mass fii, a coupling constant a^, and an 
asymptotic potential horizon Ci{T) at |r| — > cx), where the subscript i = F and U stands for the free energy or the 
internal energy respectively, 



and 




FIG. 2; The color-singlet parameters Cf, of, and jiF for the free energy _F'i(r, T) as given in Eq. 14611 . 

The parameters d, ai, and fii for the quenched lattice QCD results of Fi{r,T) and Ui{r,T) of Kaczmarek et al. 
are shown in Figs. 2 and 3, and the corresponding fits to the lattice Fi and Ui results are shown in Figs. 4 and 
5. The coupling constant apiT) is 0.44 at T = 1.05 Tc. As the temperature increases, ap decreases and saturates 
at ai? ~ 0.3 at T ^ 3Tc. The screening mass fip about 0.25 GeV at temperatures just above Tc and it increases 
to 0.8 GeV at 3Tc. The coupling constant ai;{T) is quite large at temperatures slightly above the phase transition 
temperature. At T = 1.13 Tc, ajj = 1.26. As the temperature increases, ajj decreases and saturates at au{T) ~ 0.4 at 
T ~ ATc- The screening mass fiu is small at temperatures just above Tc- As the temperature increases, the screening 
mass fijj increases to about 0.8 GeV at T ~ 4.5rc. 



The comparison in Fig. 4 and 5 shows that the free energy Fi and the internal energy Ui with the set of param- 
eters in Figs. 2 and 3, adequately describe the lattice-gauge data and can be used to calculate the eigenvalues and 
eigenfunctions of heavy quarkonia. In the local energy-density approximation, the Q-Q potential is given by 

t/W (r,T) - t/W (r ^ oo,T) = _ 4 /.(T) a.(T)e-^- + /.(T) a.(T)e-^-- ^ ^^3^ 
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FIG. 3: The color-singlet parameters Cu, ctu, and jiu for the total internal energy Ui{r,T) as given in Eq. 1471 . 
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FIG. 4: The symbols represent the quenched lattice QCD free energy, Fi{r,T)/a^^'^ , of Kaczmarek et al. |23 at selected 
distances and the curves are the fits using the screened Coulomb potential, Eq. H46^ . with parameters given in Fig. 2. Here 
o-^/^ = 425 MeV. 



VIII. CHARMONIUM IN THE QUARK-GLUON PLASMA 

In the quenched approximation, the transition temperature is Tc = 269 MeV Q . We use this value of Tc to express 
the potential in GeV units, in order to evaluate the energy levels of different heavy quarkonia. 

For a given temperature, we use the U^Q{r,T) potential given in Eq. (|43|l to calculate the charmonium energy 

levels and wave functions. In these calculations, we employ a charm quark mass rric = 1.3 and 1.5 GeV Q to provide 
an indication of the uncertainties of the eigenenergies. 

Energy levels of charmonium calculated with the u'^J.iY, T) potential are shown in Fig. 6 as a function of temper- 
ature. The J/ip and 77c states are weakly bound at temperatures slightly greater than Tc. The eigenenergies of J/ip 
and rjc are -0.045 GeV at T = 1.13 Tc for rric = 1.5 GeV and their energies increase as the temperature increases. The 
J/ip and rjc state eigenenergies are -0.0004 GeV at T = 1.65 Tc for rric = 1-5 GeV. If one extrapolates the eigenenergy 
from lower temperature points, one infers that the J/ijj and rjc spontaneous dissociation temperature is 1.52 Tc for 
■trie = 1-3 GeV, and is 1.72 Tc for rUc = 1.5 GeV, with a mean value of 1.62 Tc. There are no bound Xc "0'; ^iid rj'c 
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FIG. 5: The symbols represent the quenched lattice QCD total internal energy, Ui (r, T) /Tc, of Kaczmarek et al. [23 | at selected 
distances and the curves are the fits using the screened Coulomb potential, Eq. (1471 . with parameters given in Fig. 3. Here 
Tc = 269 MeV and a^/^ = 425 MeV. 
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FIG. 6: Energy levels of charmonium in the quark-gluon plasma as a function of temperature calculated with (a) the F\{v,T) 
and U^}l,(r,T) potentials for J/^p and 77c, and (b) the Ui(r,T) potential for J/ip,ric,Xc,i^' , and 77^. The dashed curves are 
obtained with rric — 1.3 GeV and the solid curves with rric = 1.5 GeV. Note the difference in the enegy scales in (a) and (b). 



states for temperatures above Tc. 

In the spectral function analyses of Asakawa et al. |^, ^| and Petreczky et al. |0, 0, ^|, the widths of J/ip 
begin to be broadened at ^ 1.6 Tc and the Xc states are found to be dissolved already at 1.13 Tc. The width of 
J/ 4' can be broadened by ghion dissociation g + J/ijj —f c + (5 which is presumably a possible process in the lattice 
gauge calculations in the spectral function analysis. The gluon dissociation width is however of the order of 0.05 to 
0.1 GeV as one can infer later from Section XII. In the numerical results of the spectral function analysis [Tol|. 
the width appears to be broadened by an amount significantly greater than this amount for gluon dissociation. It is 
therefore reasonable to associate the broadening of the width of a heavy quarkonium in the spectral function analysis 
with the occurrence of spontaneous dissociation, when the heavy quarkonium becomes unbound. The temperature 
at which the width of a quarkonium begins to broaden significantly can be taken as the dissociation temperature for 
spontaneous dissociation of the quarkonium in the spectral function analysis. 
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FIG. 7: J/tJj wave function calculated with the U^^{r,T) potential for different temperatures. 



In Table I we list the dissociation temperatures of different quarkonia obtained in quenched QCD. A comparison of 
the dissociation temperatures from the spectral function analysis 0,0,0,0,0 and from the u'^q{y,T) potential 
model indicates the agreement that J /ij} is bound up to about 1.6Tc and is unstable at very high temperatures. There 
is also the agreement that the Xc a-nd "0' states arc unbound in the quark-gluon plasma. 

It is instructive to compare the eigencnergics obtained by using other potential models. We calculate the heavy 
quarkonium eigenenergies with the Fi{v, T) potential as in p5ll2^l27Ll2^ by replacing C/l?! in Eq. H24|) with Fi(r, T). 
In Fig. 7(a), we show the charmonium energies calculated with the Fi(r,T) potential. One finds that J/iJj is weakly 
bound, but the dissociation temperature lies in the range 1.33-1.46 for mc ~ 1.3 — 1.5 GcV, with a mean value of 
1.40 Tc . This dissociation temperature is lower than that inferred from the spectral function analysis. The Xc and 
■0' states are unbound in this potential. 

We also calculate the heavy quarkonium eigencnergics with the total internal energy Ui{r, T) as the Q-Q potential 
as in 113, nil . The eigenenergies for charmonium states are shown in Fig. 7(b). As the total internal energy Ui 

contains a deeper potential well, the charmonium states are deeply bound. The binding energy is about 0.8 GeV at 
1.13 Tc, and the state becomes unbound at 2.50-2.71 Tc, with a mean value of 2.60 Tc- This dissociation temperature 
is much higher than dissociation temperature of about 1.6 Tc obtained from the spectral function analysis. There are 
uncertainties in the spontaneous dissociation temperatures due to the differences in the degrees of freedom assumed 
in lattice QCD calculations. For example, using the total internal energy Ui extracted from the full QCD with two 
flavors obtained by Kaczmarek et al. [TOj, Shuryak found that the dissociation temperature of J/ip is about 2.7 Tc 
|15|. In this model, the Xci "0' and 77^ states are bound at temperature slightly above Tc and they become unbound 
at 1.2 Tc- The binding of Xc states is in disagreement with the results of the dissolution of Xc states in the spectral 
function analysis. 

Our comparison of these results indicate that the model that compares best with the spectral function analysis is 
the UgQ{r, T) potential [Eq. H43|) ] obtained from a variational principle and the quark-gluon plasma equation of state. 
However, as the results from the spectral function analysis are still scanty, more results from the spectral function 
analysis are needed to test further the potential model. 



Table I. Dissociation temperatures obtained from different analyses in quenched QCD. 
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Heavy Quarkonium 


C/^^^(r,T) Potential 


Fi(r,T) Potential 


C/i(r,r) Potential 


Spectral Analysis 




1.62 Tc 


1.40 


2.60 Tc 


- 1.6 


Xc 


unbound in QGP 


unbound in QGP 


1.19Tc 


dissolved below 1 . 1 Tc 




unbound in QGP 


unbound in QGP 


1.20 Tc 






4.10Tc 


3.50 Tc 


- 5.0 




Xb 


1.18 Te 


1.10 Tc 


1.73 Tc 






1.38 Tc 


1.19 Tc 


2.28 Tc 





The solution of the Sehrodinger equation (|19|l gives both eigenenergies and eigenf unctions. We show in Fig. 8 the 
wave functions of J/i/j calculated with the U^q potential and rric = 1.5 GeV. They are normalized according to 

/"Oo poo 

/ \uis{r)fdr= / \ri:is{r)fdr = 1, (49) 
Jo Jo 

as in Eq. (4.18) of Blatt and Weisskopf [tJI. As the temperature increases, the binding of the state becomes weaker 
and the wave function extends to greater distances. The rms r of the J/tp wave function is 0.88 fm at 1.13Tc. At 
T = 1.65 Tc, which is near the temperature for spontaneous dissociation, the binding energy is 0.0004 GeV. The rms 
r of the J/ijj wave function is 5.30 fm, which is much greater than the theoretical rms r of 0.404 fm for J/^ at zero 
temperature [t^ ]. 



IX. BOTTOMIUM BOUND STATES IN THE QUARK-GLUON PLASMA 




1.0 1.5 2.0 2.5 

T/T 

c 



FIG. 8: Energy levels of bb bound states as a function of temperature, calculated with the uj^).{r,T) potential. 

WW 

One can carry out similar calculations for the energy levels and wave functions of the b-b system. We take the 
mass of the bottom quark to be 4.3 GeV, which falls within the range of 4.1 to 4.5 GeV in the PDG listing 163]. The 

energy levels of the lowest bb bound states calculated with the uj^)^{r,T) potential, are shown in Fig. 9 as a function 
of temperature. We find that at T = 1.13 T,,, the T state lies at about -0.3 GeV and the state energy increases as the 
temperature increases. The T state remains to be bound by 0.028 GeV at T = 2.5 Tc. If one extrapolates from lower 
temperatures, the dissociation temperature of T and rjb is 4.10 Tc. For this potential, the Xb, and 77^, states are 
weakly bound at temperatures slightly greater than Tc- The Xb states become unbound at 1.18X1; and the T' and rj^ 
become unbound at 1.38Tc. 
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FIG. 9: Energy levels of bh bound states as a function of temperature calculated with (a)tlie F\{v,T) potential and (b) the 
f/i(r, T) potential. 



As a comparison, we calculate the bottomium eigenenergies using other potential models. If we assume that the 
Q-Q potential is the free energy, Fi(r,T), we find that T is bound by about 0.2 GeV at l.lSTc and it becomes 
unbound at 3.50Tc as shown in Fig. 10(a). The Xb states are unbound at l.lOTc, and the T' and 77^, states are bound 
by .003 GeV at l.lSTc, which should be close to its dissociation temperature. 

If we assume that the Q-Q potential is the total internal energy, C/i(r,T), we find that the bottomium states 
become deeply bound as shown in Fig. 10(b). At 1.13 Tc, the T and r]b states are bound by about 3 GeV. The binding 
energy decreases slowly as a function of temperature. The small binding energy at T = 4.5 Tc indicates that the T 
dissociation temperature is close to and slightly greater than T ^ 5.0 Tc. Because of the small screening mass near 
Tc, the potential for temperatures near Tc is approximately a Coulomb potential but with a large coupling constant. 
Hence, Xb and T', 77^ states are nearly degenerate. They lie at -0.5 GeV for T = 1.13 Tc, and begin to be unbound at 
1.73 Tc and 2.28 Tc respectively. 

The bb bound state wave functions can be obtained by solving the Schrodinger equation H19(l . We show in Fig. 
11 the T radial wave functions as a function of temperature. The wave function extends to greater distances as the 
temperature increases. The root-mean-square radius rims is 0.25 fm at T/Tc = 1.13 and it increases as the temperature 
increases. At T = 2.6 Tc, which lies very close to the temperature for spontaneous dissociation, the rms r value is 0.59 
fm. which is substantially greater than the rms r of 0.25 fm at 1.13 Tc. 

Our comparisons of the charmonium eigenenergies with those from the spectral function analysis indicate that the 
results obtained by using the potential give results which agree best with the spectral function analysis. As 

different models give different predications on the bottomium dissociation temperatures, further tests of the models 
can be carried out by calculating the dissociation temperature of bottomium states using the spectral function analysis. 

X. ANTISCREENING BY DECONFINED GLUONS 

As one crosses the phase transition temperature Tc from below, quarks and gluons becomes deconfined. The Debye 
screening due to the interaction between the heavy quark and antiquark with medium particles is considered to be the 
dominant effect when a heavy quarkonium is placed in a quark-gluon plasma 0. It is argued that Debye screening 
leads to a decrease in the attractive interaction between the heavy quark and antiquark and results in the spontaneous 
dissociation of the heavy quarkonium in the quark-gluon plasma. The suppression of J/ijj production was suggested 
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FIG. 10: T wave function rilj{r) as a function of temperature. 



as a signature for the quark-gluon plasma . 

In perturbative QCD, there is a relation between the screening mass /i, coupling constant g = v47rasj temper- 
ature T in a quark-gluon plasma given by [7^ O IzE lzl> Iz3 

= .g2r2(iVe + Af//2)/3. (50) 

For the quenched approximation with N f ~ 0, perturbative QCD gives /i ~ gT. A comparison of the effective coupling 
constants ap and au of either Fi or Ui with the screening mass fip and yU(7 at temperatures near Tc indicates that 
the screening mass is much smaller than the estimates in perturbative QCD limit, indicating that the large-distance 
behavior cannot be qualitatively described by perturbative QCD. The perturbative QCD limit for the large-distance 
behavior can be reached only at temperatures above 6 Tc [t^ l79l|. In fact, Kaczmarek et at. |79l | pointed out that 
for temperatures close to the QCD phase should be more appropriately described in terms of remnants of the 
confinement part of QCD rather than a strongly coupled Coloumbic force |14l Il5l|. To understand this "remnant of 
QCD confinement above Tc" , it is of interest to examine effects of antiscreening above by deconfined gluons. 

In a related topic, Svetitsky, Yaffe, DeTar and DeCrand previously found that large space-like Wilson loops in the 
quark-gluon plasma have an area law behavior 0, 0, Is^ Isp . Is^ l . As the area law of the spatial Polyakov loop does not 
change drastically across the phase transition, DeTarj^jU, Hansson, Lee, and Zahed and Simonov |EIEI3 argued 
that low-lying mesons including J /ip may remain relatively narrow states due to the attractive interaction between 
the quark and the antiquark and the suppression of J /ip is not a signature of the deconfinement phase transition |3| ■ 
It is of interest to examine the consequence of the spatial area law to see whether it will lead to antiscreening between 
a heavy quark and antiquark. 

The mechanism of antiscreening by virtual gluons at T = is well known (see for example Peskin, Schroeder, 
Gottfried, and Weisskopf 113, 1^)- We would like to follow similar arguments to study the mechanism of antiscreening 
by deconfined gluons above the phase transition temperature Tc- We consider the color electric field generated by a 
static color source p°^j''(r) = 5{y)5'^^ with a unit color charge of index A placed at the origin, in the presence of an 
external gauge field (Fig. 12a). We fix the gauge to be the Coulomb gauge and represent the deconfined gluons in 
terms of an external transverse gluon field A^'^{y) where the first index 5 = 1, 8 is the color index and the second 
index i = 1, 2, 3 is the spatial coordinate index. 

The color electric field i?"' (r) generated by the source is determined by the Gauss law 

^^E''\v) = gS{r)S''^ + gf'^^ A^\v)E"{r) (51) 

(Eq. (16.139) of Peskin and Schroeder [sJl). Here repeated indices are summed over and the first index of i?°*(r) is 
the color index and the second index is the spatial coordinate index. Because of the non-linear nature of the second 
term which arises from the non-Abelian nature of QCD, the external color source 5{y)5°-^ and the external gauge field 

A°*(r) induce a color source p'^i^i'^)^ which in turn induces an additional color source p^^i^). How do these induced 
color charges depend on the external gauge field A^^{y)1 

We consider an expansion of the source in terms of the external source and the induced sources, in powers of the 
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coupling constant 

d.E^\r) = .9<5(x)5°i +5Pr„d^W +5Pr„d^(r) + ... (52) 
In the Coulomb gauge, the color field i?'^*'^^^ (r), as arising only from the external static source S{r)6°'^, is 

(r) =g,5°^^. (53) 



From the non-linear term in Eq. (|51|l . the color charge density induced at r by the external gauge field A'^*(r) and 
the electric field E'^^^^\r) of the external color source is 

pfndV) = r^^AP\r)E^^^^\r) = .g/-/37^/3'(r)^. (54) 

An induced color-charge element p°JP(r)Ar at r will generate a field E°-'^'^'^\r' , r) at r' and this field is pointing along 
the direction of r' — r, 

£;'^^(2)(r',r)=g[pr«(r)Ar]iJ^. (55) 

From the non-linear term in Eq. (|51|1 . the color charge density element p'^^'^\r' , r) Ar , that is induced at r' by the 
external gauge field A'"(r') and the electric field (r', r) at r', is therefore 



t:^- (56) 



As the external source has the color index A, we would like to study the induced color charge of index A to see whether 
the induced color charges lead to screening or antiscreening. From Eq. (|54|) . we have 

P-!d'^(r) = 0, (57) 

as f^^'^ = on account of the antisymmetric property of /. For the next-order induced color charge density element 
(r', r)Ar, we can write out explicitly the summations of color and spatial indices of Eq. H56|l . 

b,cJ3=l ij = l ' ' ' ' 

We note that 

8 
c=l 

where i^(A, b) is a non-negative quantity defined by 

F(A,&)=^(/^^^)^ (60) 

c=l 

which can be easily evaluated. For example, F{l,b) is {0, 1, 1, 1/4, 1/4, 1/4, 1/4, 0} for b — 1,2,. ..8. In terms of 
F{X, b), the induced charge density is 

Note that F{X, A) = 0. The contribution of the external gauge field to the sum over b comes only from those color 
components of A*" that are transverse to the color axis A of the external point source. 

The product A^^{r')A''^ {r) in Eq. (|61|l involves spatial gauge fields at different spatial locations. Previously, Svetit- 
sky, Yaffe, DeTar and DeGrand found that large space-like Wilson loops in the quark-gluon plasma have an area 
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FIG. 11: (a) The coordinate system used in the evaluation of the induced color charge density pj^^ (r',r) at r'. (&) A heavy 
quark of color source at O induces color charges of the same sign in the direction forward of r and color charges of the opposite 
sign in the direction backward of r. 



law behavior. Such a behavior indicates that the spatial gauge fields A^^ (r) at different spatial locations are corre- 
lated 0, H, 1^ Isl, 01 . As shown in Appendix A, if the gauge fields A^"^ at different field points r> and r< with 
|A''*(r>)| > \A^''{v<)\ are correlated as 

A^\y<) = ^''*(r>)e-l"">"'"<l/«, (62) 

then the integral of A^^(t) along a space-like loop § A^'^dx^ obeys an area law, when the linear dimensions of the loop 
are small compared with the correlation length. Consequently, within the vicinity of the the correlation length ^, we 
can express the thermal average of the A^'^{y')A^^ [r] in terms of the relative coordinate and the correlation length ^. 

{A^\y')A^^{y)) - A''*(r)A''^(r)e-l'"'-'"l/«. (63) 

where we have considered the case when |A''*(r')| < |A''*(r)|. (The case of |A''*(r')| > |A'''(r)| can be treated in 
an analogous way). The quenched QCD calculations in 5?7(3) is in the universal class of three-dimensional Z{i) 
symmetr y. It has a first ord er p hase transition and and it possesses a large but finite correlation length at 

II II 111 M in m m m m ii . 

To study the sign of the induced color charges, we choose a spatial coordinate system with the z-axis along r as 
shown in Fig. (4a). In this coordinate system, we label the angular coordinates of r' — r and A''(r) by {0',(j>') and 
{db,4>b), respectively. We shall study the induced charge at r' within the vicinity of the correlation length ^ of r. 
Then, after taking the thermal average, the induced color charge density is then 

6=1 ' ' 

X cos 0b [cos 6*6 cos 0' + sin6'b sin(?' cos((/>(, — </)')]. (64) 



If we average over the azimuthal angle 0', the second term in the square bracket drops out and we have 

^ |A''(r)|2 cos^eifc e-l""'-''l/«cos6i' 

j,2f^j.2 _ j.|2 ^ 2r\r' — r| cos 6') ' 



\(2), , ^^ 9V^^/, A" r)r cos'^h e^i' "'I'^cosy 

h=l 

One readily observes that the induced color charge density p^^^ {r' ,r) at r' is negative in the forward hemisphere in 
the direction forward of r (with 7r/2 > 6*' > 0). It changes to positive in the backward hemisphere, in the direction 
backward of r (tt > 0' > 7r/2). In the region of r' within the vicinity of the correlation length from r, the induced charge 
surrounding r is a color-dipole type density distribution with the color charge of the same sign at distances closer to 
the color source and of the opposite sign at distances farther to the color source (Fig. 46). This is the antiscreening 
behavior due to the presence of external gauge field A^'^{r) at r. The magnitude of the induced color charges will 
increase with an increases in the correlation length ^ and the magnitude of the gluon field. The antiscreening effects 
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will enhance the attractive interaction between the heavy quark and antiquark and will reduce the screening mass 
from the Debye screening predictions. 

The antiscreening effect arises due to the non-linear properties of the non-Abelian gauge field while the effects of 
Debye screening arises from the interaction between the quark and antiquark with gluons. Both effects are present 
and the antiscreening effects due to deconfined gluons will act to counterbalance the effects of Debye screening. At the 
onset of the phase transition, the correlation length is large [111,!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! and deconfined 
gluons are present, there can be "remnants of the confinement part of QCD" at temperatures slightly above Tc, as 
pointed out by Kaczmarek et al. [t^- At a much higher temperature, a greater thermal fluctuation leads to a smaller 
correlation length, reduces the effects of antiscreening, and Debye screening dominates. 



XI. DISSOCIATION OF J/iP IN COLLISION WITH GLUONS 

It is not necessary to reach the spontaneous dissociation temperature with zero binding energy for a quarkonium to 
dissociate. In a quark-gluon plasma, gluons and quarks can collide with a color-singlet heavy quarkonium to lead to 
the dissociate of the heavy quarkonium. Dissociation by the absorption of a single gluon is allowed as the color-octet 
final state of a free quark and a free antiquark can propagate in the color medium, in contrast to the situation below 
Tc in which the quark and the antiquark are confined. We shall consider dissociation of heavy quarkonium by gluons 
in the present work. The collision of a heavy quarkonium with light quarks can also lead to the dissociation of the 
heavy quarkonium, but through higher-order processes. They can be considered in future refined treatment of the 
dissociation process. 

Previous treatment of the dissociation of heavy quarkonium by the absorption of a gluon was carried out by Peskin 
and Bhanot [s^lo^- They use the operator product expansion and the large Nc limit. They sum over a large set of 
diagrams and show that to obtain gauge invariant results, they need to sum over diagrams in which the external gluon 
is coupled to the gluon that is exchanged between the heavy quark and the heavy antiquark. They use hydrogen wave 
function and hydrogen states to evaluate the transition matrix elements. Their expression for the dissociation cross 
section of a{g + {QQ)is Q + Q) is 

+ (00)« ^ + 0, 4, (f ) ' (^) -L A^^, (ee, 

where E is the non-relativistic kinetic energy of the dissociated Q and Q in the center-of-mass system. In this short- 
distance approach of Peskin and Bhanot, the quark and the antiquark form a color-dipole pair and the gluons couple 
to the Wilson loop (the quarkonium) via simple dipole interactions. The dissociation cross section of Eq. H66|) is, in 
fact, the dissociation of the quarkonium through the absorption of an El gluon radiation. 

Peskin and Bhanot 's analytical result for the dissociation cross section has been applied to many calculations 
[23L I23 . l9ll|. In heavy quarkonia of interest, the radial dependence of the quark-antiquark potential often differ from 
the Coulomb potential. The calculation of the dissociation cross section requires a new formulation which can be 
best described by the potential model introduced previously |2^, , following the results of Blatt and Weisskopf 
[tH obtained for the photo-disintegration of a deuteron. The dissociation process is schematically illustrated in Fig. 
12. An initial bound {QQ)is state with a binding energy B in the color-singlet potential Vi{r) absorbs a gluon of 
energy hujg, and is excited to the color-octet final state {Q + Q)ip with a kinetic energy E above the rest mass of 
ruQ + ruQ. The interaction V3{r) between Q and Q in the color-octet state will be different from the interaction 
Vi(r) in the color-singlet state, as shown in Fig. 12. At low energies, the dominant dissociation cross section is the 
El color-electric dipole transition for which the final state of Q -I- Q will be in the IP state in the continuum. The 
dissociation cross section (j{g + J/ih c + c) for such a color El transition can be obtained from the analogous result 



in QED lilllJl, and the result is H 



where 



^<fis(^giuon) = 4 X -agQ(/c2 + 72)fc-i J2, (67) 



7^ ^ 2nB, fc2 = 2nE, (68) 
I ^ I uip{r) r uis{r)dr, (69) 







«.|(8c|y|l)|^ = a. X 1, (70) 
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FIG. 12: The quarkonium dissociation process in the potential model. 



and tts is the gluon-(heavy quark) coupling leading to the dissociation of the heavy quarkonium. Here, we use the 
same notation and normalization as in Blatt and Weisskopf. The bound state wave function uis has been normalized 
according to Eq. H49|) as in Eq. (XII. 4. 18) of Blatt and Weisskopf [tJ, and the continuum wave function uip is 
normalized according to 



Uip{r)\r 



krji{kr) 



sin(fcr) 
kr 



cos(fcr), 



(71) 



as in Eq. (XII. 4. 32) of Blatt and Weisskopf [71|. The result from the potential model agrees with the analytical 
results of Bhanot and Peskin for the case they considered (hydrogen wave function, large Nc limit,...) |2^. Such 
an agreement was further confirmed by numerical calculations according to Eq. H67|) using h ydr ogen wave function 
for Mis and plane wave continuum wave function for uip, as assumed by Peskin and Bhanot |89|, |90j. The potential 
model has the practical advantage that it can be used for a Q-Q system with a general potential. 



XII. J/ip COLLISIONAL DISSOCIATION RATE AND DISSOCIATION WIDTH 

We have represented the color-singlet potential between a heavy quark and antiquark by a screened Coulomb 

potential and have obtained the J/^j wave function. To study the gluon dissociation of J/iJj, we need the color-octet 

(s) 

[)()! cut iai V8(r, T) = U^q{t,T) experienced by the Q and Q in the final state. We shall assume the generalization 
that the color-dependence of the potential Eq. (|46|) is simply obtained by modifying the color factor from —4/3 for 
the color-singlet state to 1/6 for the color-octet state; 

U^ir^T)-U^ir^o..T) = c /^"- ^"'^ ^ ^^"^ ^'^^ \ (72) 

C = /"^/^ (i = 1, color -singlet) 
^ y 1/6 (i = 2, color- octet). ^ ' 

We also need the gluon-quark coupling constant a,, in Eq. H70() to evaluate the dissociation cross section. Wc shall 
consider the screened Coulomb potential (|67|) obtained in the lattice gauge theory as arising from the exchange of 
a virtual nonperturbative gluon and assume that the coupling of gluon to the heavy quark leading to quarkonium 
bound states, fpoip + fuctu, is the same coupling leading to the dissociation of the quarkonium. The J/tp dissociation 
cross section can then be calculated using Eq. H67|l . The results of the J/-0 dissociation cross section for different 
temperatures are shown in Fig. 13. The cross section increases up to a maximum value and decreases as the gluon 
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energy increases. As the temperature decreases the maximum height of the dissociation cross section increases, but 
the width of the cross section decreases. We shall limit our attention to the dissociation of J /ip by the absorption of 
an El radiation in the present analysis. When the dissociation threshold decreases, higher multipole dissociation may 
become important. It will be of interest to study dissociation arising from gluon radiation of higher multipolarity in 
future work. 




gluon 



0.3 0.4 

(GeV) 



FIG. 13: J/tp dissociation cross section as a function of gluon energy at various quark-gluon plasma temperatures. 



If we place a J/tp in a quark-gluon plasma at a temperature T, the average El dissociation cross section is 



.Bix _ 9g ^Elt^^_P dp 

e 



dis / - ^ / ^dis (P) ^p/T „ ^ / "9 ' (74) 



where gg — 16 is the gluon degeneracy and Ug is the gluon density, 

_ 9g P dp 



271-2 / gp/T _ I ■ 



(75) 



Using the energy dependence of the dissociation cross section as given in Fig. 13, we can calculate the average 
dissociation cross section (crfjg). From (crfjg), we obtain the rate of J ftp dissociation (by El transition) given by 

-ngia^-J. (76) 

This dissociation rate leads to the coUisional dissociation width Fei due to the absorption of El gluon radiation, 

rEi= ngia^). (77) 

We show in Fig. 14 the temperature dependence of (crfjs) and Fgi. One observes that the average cross section is 
in the range of 0.2-0.8 mb. The coUisional dissociation width due to El gluon absorption is of the order of 0.05-0.11 
GeV, and the mean life of J/^ in the quark-gluon plasma due to the absorption of gluons to the IP state is therefore 
of order 2-4 fm/c. 



XIII. J/i: PRODUCTION BY THE COLLISION OF c AND c IN QGP 



In high-energy nuclear collisions, elementary nucleon-nucleon collisions lead to the production of open heavy quark 
mesons. Although the probability for such a production is small for a single nucleon-nucleon collision, there are many 
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FIG. 14: (a) Thermally averaged dissociation cross section as a function of temperature, (b) J/t/) coUisional dissociation width 
r_Bi as a function of temperature. 



nuclcon-nuclcon collisions in a central nucleus-nucleus collision. The large number of binary nueleon-nucleon collisions 
can produce many pairs of open charm mesons, and these open charm mesons can recombine to produce J l'4>- It is of 
interest to estimate the elementary reaction cross section for c + c — *■ J/?/; + g and obtain the rate of J/?/; production 
in a nucleus-nucleus collision. 

The reaction c + c ^ J/V' + 5 is just the inverse of g + J/V' ^ c + c. Their cross sections are therefore related by 

^^'(c + c J/^ -f g) = j/v,, + g ^ c + c), (78) 

where is the momentum of one of the particles in the J /i}}+g system, and is the momentum of one of the particles 
in the c-|-c system, both measured in the center-of-mass frame. With Eq. (|78|) and the results of (T^^{J/ip + g c + c) 
in Fig. 13, the production cross section a^^{c + c ^ J /ip + g) can be calculated. The cross section as a function of 
the kinetic energy of c and c in the center-of-mass system are shown in Fig. 15. One observes that the cross section 
peaks at low kinetic energies near the threshold, and the magnitude of the cross section decreases with temperature. 
The maximum cross section at T/Tc ~ 1.13 is of order 0.7 mb. 

The rate of J/ip production can be obtained when the momentum distribution ,f{y,Pt) of the produced c and c 
is known. For simplicity, we consider charm quarks and antiquarks to be contained in a spatial volume V with a 
uniform distribution. The probability of producing a J /ip in the volume V per unit time by the collision of charm 
quark and antiquark is cr(c + c ^ J /ij) + g)vi2/V, where 1112 is the relative velocity between c and c. The number of 
J/ip produced per unit time from collision of c and c is therefore 



j dyidpitdy2dp2tf{yi,Pit)fiy2,P2t)(^(c + c-> J/'iij + g)vr2/V. (79) 



dt 

When we neglect initial- and final-state interactions, the momentum distribution of charm is given by 

/(y,pJ=iVbi„^^. (80) 

Here /(y, p^) is normalized to the total number of charm quarks Nc and antiquarks Ng produced in the nucleus-nucleus 
collision, Nc = Nc = J {dp / E) f {y , p^) = A^bin^ce': ^bin is the number of binary nueleon-nucleon coUisions, and Nc^ is 
the number of cc pair produced in a single nueleon-nucleon collision. From the charm production data in d-Au collisions 
and PYTHIA calculations, Tai et al. [s^i Adams et al. and the STAR Collaboration inferred that at ^/s = 200 
GeV the charm production cross section per nueleon-nucleon collision is cr^g |star = 1-18 ± 0.21(stat) ± 0.39(sys) 
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FIG. 15: Cross section for the production of J/ip by the coUision of c and c 



mb [23 E^nd cTceIstar = 1-4 ± 0.2 ± 0.4 mb If one uses the transverse momentum distribution measured and 

parametrized by Tai et at. j9l^ and assumes a Gaussian rapidity distribution, the charm momentum distribution of 
Adams et al. |9J| per nuclcon-nucleon coUision can be represented by 



EdN^ 
dp 



(PYTHIA(STAR)) 



STAR 

aindp 



2 /o 2 



A 



(81) 



where Gin = 42 mb is the nuclcon-nucleon inelastic cross section, A = 4.4 x 10"'^ (GcV^^), Uy = 1.84, pto = 3.5 GcV, 
and n = 8.3. The number of cc produced per nucleon-nucleon collision is A^^^^'^Istar = 1-4 mb/42 mb=0.033±0.0107. 

The PHENIX Collaboration obtained cTcflpHENix = 622±57(stat)±160(sys) fih [93 for the open charm production 
cross section per nucleon-nucleon collision at -y/s = 200 GeV, and A'^g^jpHENix = 0.0148±0.004. With this total charm 
production cross section, the theoretical results from the PYTHIA calculations of the PHENIX Collaboration can be 
parametrized as j95|, l96l] 



EdN^^ e-y^^^^y 
^"''--.(PYTHIACPHENIX)) = A'- 



dp 



(l+Pt/Pt 



(82) 



where A' = 6.48 x 10"'^ (GeV'^), a'^ = 1.85, pJo = 5.06 GcV, and n' = 7.0. 

We shall focus attention to central Au-Au collisions within the most inelastic 10% of the reaction cross section. The 
average number of binary collision A^bin for these central Au-Au collisions at RHIC is A^bin = 833. For these central 
nucleus-nucleus collisions at ^/s = 200 GeV, the average number of c and c produced is Nc = Ns = 27.8 ± 8.9 if we 
use the cross section of the STAR Collaboration |93|, and = Nc = 12.34 ± 3.4 if we use the cross section of the 
PHENIX Collaboration The rate of J/^p production can then be obtained from Eqs. (|79|I - H82|I by carrying out 
the six-fold integration. 

In Fig. 16 we show the quantity VdNj/^/dt as a function of temperature for the most inelastic (10%) central 
Au-Au collisions at = 200 GeV. The estimate of the rate of J /ip production, using the momentum distribution of 
Eq. H81|) from the PYTHIA calculations of the STAR Collaboration [9J|, is greater than the corresponding estimate, 
using the momentum distribution of Eq. from the PYTHIA calculations of the PHENIX Collaboration , 
by a factor of about 10. This factor is larger than the factor of 2.25 of the nuclcon-nucleon cc production cross section 
of the STAR Collaboration, relative to the corresponding cross section of the PHENIX Collaboration at -^s ~ 200 
GeV. The large difference of these two factors arises because the charm momentum distribution from the PYTHIA 
calculations of the STAR Collaboration has a greater magnitude at low py, and c and c recombine more readily at 
low relative energies, as indicated in Fig. 16. 
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FIG. 16: The rate of J/i}i production as a function of temperature, for the most inelastic (10%) central Au-Au collision at 
y/s = 200 GeV. Curve A is based on the charm momentum distribution of Eq. I|81|l using the PYTHIA calculations of the STAR 
Collaboration loj . and Curve B is based on the charm momentum distribution of Eq. 18211 using the PYTHIA calculations of 
the PHENIX Collaboration 

We can illustrate the magnitude of the rate of J /ij) production by considering a Au-Au central collision with a 
transverse area of 7r(7 fm)^ and a longitudinal initial time (and longitudinal length) of 1 fm. The initial volume 
containing the charm quarks and antiquarks is about 150 fm^. If the initial temperature is lATc, the initial rate of 
J/i/j production will be about 3 x 10^^ to 3 x 10"'^ fm/c. As the volume expands, the temperature decreases. The 
results of Fig. 17 can be use to provide an estimate of the rate of J/V" production. 

XIV. DISCUSSIONS AND CONCLUSIONS 

We use the color-singlet free energy Fi and internal energy Ui obtained by Kaczmarek et al. in quenched QCD 
to study the energy levels of charmonium and bottomium above the phase transition temperature. From a variational 
principle in a schematic model, we find the the Q-Q potential involves only the Q-Q internal energy i/lU, which can be 
obtained from the total t/i by subtracting the gluon internal energy contributions. We carry out this subtraction using 
the local energy-density approximation in which the gluon energy density can be related to the local gluon pressure 
by the quark-gluon plasma equation of state. We find that the Q-Q potential is C/gg = 3Fi/(3 + a) + aUi/{3 + a) 

where a ~ 3p/e is given by the quark-gluon plasma equation of state. Such a f/lU potential leads to weakly bound 

J/ip and rjc at temperatures above the phase transition temperature and they become unbound at 1.62Tc. The Xc, 
r]'^ and tp' states are found to be unbound in the quark-gluon plasma. In this potential model, T, rjb, T', 77^ and x& 
are bound at temperatures above Tc and T and rjb dissociates spontaneously at 4.1T'c, Xb at l.lSTc and T', and 77^ at 
1.38Tc. 

The results from the U^q potential need to be tested against results from spectral function analysis. For complete- 
ness, we have also calculated heavy quarkonium binding energies using the free energy Fi |2^ |2^ |2^ and the 
internal energy Ui as the potential |3.l29ll30j. 

The comparison shows that different models give very different heavy quarkonium binding energies. The potential 
that agrees best with results obtained from spectral function analysis is the U^q potential deduced in the present 
analysis. The agreement with spectral function analysis and the theoretical foundations presented here provide support 
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for the use of UXA as the proper Q-Q potential in heavy quarkonium studies. Conversely, the agreement also lends 
support to the quantitative features concernin g th e stability of heavy quarkonia in the spectral function analyses of 
Asakawa et al. and Petreczky et al. [TllTTllT3|. 

The spectral function analysis for the bottomium states has not yet been carried out. As the predications for the 
dissociation temperatures for bottomium states are quite different from different potential models, it will be of great 
interest to calculate the bottomium dissociation temperatures in lattice gauge spectral function analysis so as to test 
the potential models further. 

In a nucleus-nucleus collision, charm quarks and antiquarks arc produced in hard-scattering processes in nucleon- 
nucleon collisions. During the time of a central nuclear reaction, these heavy quarks and antiquarks will be present 
in the quark-gluon plasma and can interact to form J/ip- We have calculated the cross section for J/ip production by 
the collision of charm quark and antiquark. The cross section is energy dependent, and the maximum cross section 
increases as the temperature decreases. The production cross section can be used to study the rate of J/ip production 
in nucleus-nucleus collisions. 

We have carried out the investigation using the quenched QCD. It will be of interest to carry out similar investiga- 
tions using unquenchcd QCD. Results of the full QCD in two flavors [t^I and in three flavors ^| have been obtained 
recently, and an investigation on J/ip dissociation temperatures in QCD with two flavors has been initiated |l5l |. A 
thorough study of how the dynamical quarks will affect the stability, the dissociation, and the inverse production of 
heavy quarkonium will be of great interest. 

It is necessary to emphasize that the present Q-Q potential C/IU of Eq. (|43|) extracted from Fi and Ui has been 
obtained in the local energy-density approximation. It would be of great interest in future lattice gauge calculations 
to evaluate the C/l|'^(r) directly to check the validity of the local energy-density approximation. 

The color-singlet correlator of the Polyakov lines in Eq. ^ is not gauge invariant. It has been suggested that one can 
dress the Polyakov lines to make a gauge invariant definition of the color singlet potential. The dressing of the source 
may be viewed as a gauge transformation and is equivalent to the choice of a certain gauge |98j with the requirement 
that the gauge-fixed Polyakov loop correlation function in the singlet channel falls off with gauge- invariant eigenvalues 
of the Hamiltonian. This requirement may be satisfied for the Coulomb gauge and other time-local gauges. Recent 
calculations by Belavin et al. show however that the color singlet potential depends on the choice of the gauges 
even among these time-local gauges. It has been found that at finite temperatures all channels receive contributions 
only from the color-singlet channel. The extraction of the color-singlet potential from the "color-singlet" Pol yako v 
correlator of Eq. |^ may include additional r and/or T dependence which is not shared by the physical states |lOCll |. 
Clearly, much work remains to be carried out to clarify the proper color-singlet potential in lattice gauge calculations. 



The authors would like to thank Drs. H. Crater, P. Petreczky, M. Gyulassy, Su Houng Lee, Keh-Fei Liu, T. Barnes, 
S. Ohta, V. Cianciolo, D. Silvermyr, Huan Huang, and Zhangbu Xu for helpful discussions and communications. 
This research was supported in part by the Division of Nuclear Physics, U.S. Department of Energy, under Contract 
No. DE-AC05-00OR22725, managed by UT-Battelle, LLC and by the National Science Foundation under contract 
NSF-Phy-0244786 at the University of Tennessee. 

APPENDIX A: SPACE-LIKE AREA LAW AND THE CORRELATION OF GAUGE FIELDS 

We focus our attention on the x ^ y plane so that the z coordinate can be omitted and consider a loop integral 
Aidx^ along the loop L defined by (0, 0) {Lx, 0) (Lx, Ly) (0, Ly) (0, 0). The integral around this loop of 
area L^Ly is given approximately by 
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We can write the right hand side in the form 




(A2) 



where 



A,{L.j2,0)-A,{L,/2,Ly) 



Ay{Lx,Ly/2)~Ay{0,Ly/2) 



(A3) 



c 
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If the gauge fields A at different field points are correlated by a correlation length ^ such that for two points r> and 
r< where |Aj(r>)| > |Ai(r<)| and 

A,(r<) = A,(r>)e-l"^>"'^<l/«, (A4) 

then we have 



A^{Lj2,Q)-A^{Lj2,Ly) _ \ {l - e-^y/^)A.,{if,0)l Ly if 0) > L,), 

Ly \{e'^y/^ ~l)A,{if,Ly)/Ly if A,{if,Q)<A,{if,Ly) 



(A5) 



The second term in Eq. (|A3() can be similarly evaluated. In the case of a correlation length ^ that is large compared 
with the loop lengths and Ly, the quantity c in Eq. (|A2|) can be evaluated and we obtain the area law 

£ Adx' = j{A>- A>)L,Ly, (A6) 



where 



^> ^ U,(i,/2,0) if A,iLj2,0)>A,iL,/2,Ly), 
" \-A4L,,/2,Ly) if A4L.,/2,0) < A4L,/2,Ly), ^ ' 

and 

^> ^ fA,(0,L,/2) if AyiO,Ly/2)> AyiL,,Ly/2), 

y yAy{0,Ly) if Ay{0,Ly/2)<Ay{L,,Ly/2). 

Eq. I|A6|I shows that if space-like gauge fields at different points are correlated by Eq. (|A4|) with a large correlation 
length, the integral of the gauge fields along a space-like Polyakov loop will satisfy an area law. 
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